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Abstract 

The model for hybrid excitations of the QCD string with quarks 
is presented starting from the perturbation theory in the nonpertur- 
bative background. The propagation of a system containing q'^-pair 
and gluon is considered. The simplified version of the Hamiltonian, 
including both long-range nonperturbative interaction and Coulomb 
force, is derived. The masses of the lowest qqg hybrids are evaluated, 
and numerical results for the spectra are listed. 



1 Introduction 

One of the most important features of the QCD Lagrangian is the presence of 
gluonic degrees of freedom which should exhibit themselves at the constituent 
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level, namely in the form of glueballs and hybrids, i.e. non-gg exotics. Un- 
fortunately, the current experimental situation is too complicated to give us 
unambigious proof that such states exist; some candidates for non-qq exotics 
appear and disappear from time to time, changing their masses and quan- 
tum numbers (for the up-to-date review see [1]). Nevertheless, there is no 
doubts that standard qq nonets arc ovcrpopTilatcd, but the question about 
the nature of the "extra" states is still open. On the other hand, there is no 
substantial progress in the description of the strong coupling nonperturba- 
tive regime in QCD. There are some QCD - inspired theoretical approaches, 
but none of them are able to provide reliable enough predictions for masses 
and decay rates of hybrid mesons. 

The QCD sum rules estimations for exotic hybrids are rather unstable: 
first results [2] predicted for the 1~+ hght hybrid the mass 1.2-1.7 GeV, while 
more recent calculations give 2.1 GeV [3] and 2.5 GeV [4]. 

In the bag model the gluons are automatically transverse, and the lowest 
electric gluon (with = l"*") in the spherical cavity is much lighter than 
the lowest magnetic one (with = 1^). The lowest hybrids with light 
quarks have the mass about 1.5 GeV [5,6]. The masses of hybrids with 
heavy quarks were estimated in [7] taking into account the bag deformation, 
with the results 3.9 GeV for cc hybrid and 10.5 GeV for bb one. 

Constituent gluon model was introduced in [8,9]. In this model the linear 
potential is introduced ad hoc, in analogy with the charmonium system. 
It was proposed also that the gluon orbital momentum is diagonal. As a 
consequence, in such model the lowest states with the S'-wave gluon have 
non-exotic quantum numbers. The mass predictions of constituent model 
are 1.3-1.8 GeV for the lowest non-exotic hybrids. 

Flux-tube model [10] predicts degenerate (up to spin corrections) light- 
est hybrid states at 1.8-1.9 GeV. In this model phonon-type excitations of 
the string connecting quark and antiquark are interpreted as hybrids. First 
calculations [10] assumed small oscillation approximation, and recent results 
[11] demonstrate that this approximation might be inadequate. Improved 
version [11] is given in the framework of "one-bead" fiux-tube model. 

Here we present the studies of the qqg system in the framework of Vacuum 
Background Correlator method [12]. The main assumption is that nonper- 
turbative background fields {B^} exist in the QCD vacuum, which ensure 
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the area law asymptotic for the Wilson loop along the closed contour C, 

< W{C) >B^ Ncexp{-aS), 

with S — Sjnin being the minimal surface inside the contour C. It was con- 
firmed by cluster expansion method [12] that if one assumes the existence 
of finite correlation length for the background, the asymptotical behaviour 
of the Wilson loop average is compatible with the area law. The deviation 
from the area law at large distances with S = Smin are caused by the pertur- 
bations over the background. The "minimal" area law gives the string-type 
interaction in the qq system, while the perturbative fields are responsible for 
the string vibrations [13]. 

2 Green function for the qqg system in the 
Vacuum Background Correlator method 

The constituent gluon in the Vacuum Background Correlator method is in- 
troduced as a gluon propagating in the nonperturbative background field 
[14,15]. Following [15], we split the gluonic field into the background field 
Bn and the perturbation over the background. 

A^^B^ + (1) 
We ascribe the inhomogeneous part of gauge transformation to the field 

B^^U+iB^+'-d^)U, a^^U+a^U, (2) 

so that the states involving the field may be formed in the gauge-invariant 
manner. One-gluon hybrid is represented as 

Xq, Xg) = i)a{Xq)^''p{Xq, X g)a'^^{x g)^1{x g, Xq)ijj\xq), (3) 

where a...S are the colour indices in the fundamental representation, = 
Oa(Aa)!^, and parallel transporters $ contain only background field: 

$^(x,y) = (Pexp fB.dzXp (4) 
Jy 
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The Green function for the qqg system is obtained by averaging the prod- 
uct ^in^aut over the background field configurations: 

G{XqXqXg, yqyqyg) =< *m(yq, Vq, yg)'^l^i{XqXqXg) >B ■ (5) 

The dynamics of the field is defined, in accordance with the decompo- 
sition (1), by expanding the QCD Lagrangian up to the second order in the 
fields (in the Euclidean space-time) 

L(a) = -^{F^,{Bjf + a,D^{B)F^,{B)+ (6) 

+^a^{DxDxS^^ - Df,D^ - gF^^{B))a^, 

Dl^^dxS^'^ + gr'^Bl 
with the background gauge fixing term 

G« = d,a; + gr'%a; (7) 

In what follows we skip the issue of ghosts. The linear in the fields 
part of the Lagrangian (6) disappears if the field B^ satisfies the classical 
equation of motion D^Ff^,^ — 0. To be on the safe side one is to assume that 
the background is the classical one, or at least, that the transition vertex 
generated by this term is small. 

The Green function for the field propagating in the given background 
B^ may be identified in the background gauge as 

G^J = D\, - D,D^ - gF^, + -^D^D, = M^, + -^D^D,. (8) 

If the classical equations of motion are satisfied, then one has M^i^D^ — 0, 
and the Green function (8) may be rewritten as 

G,. = {S^x + (e - l)D,^Dx){D' - 2gF)-xl. (9) 

The choice ^ = corresponds to the Landau gauge, in which the Green 
function (9) contains explicitly the projector P^x onto transverse states: 

P,x = 5,x - D,{^)Dx. (10) 
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To define the effective action for the qqg system we use the Feynman- 
Schwinger representation [16]. To do it in proper way one should talce into 
account spin degrees of freedom of quarks and gluon. Here we omit the spin 
dependence, reducing the problem to the scalar one. This simplified version 
of the model corresponds to the neglecting of colour magnetic interaction (the 
term proportional to the gF in (9)), and omitting the projector. Similarly, 
spin dependence in the quark Green function is also omitted, and we assume 

G,^{D'-ml)-\ (11) 

As the result, the Feynman-Schwinger representation for the hybrid Green 
function takes the form 

G{XqXqXg, yqVqyg) = (12) 

poo poo poo _ p 

ds ds dS j DzDzDZexpi-JCg - -]Cg) <W >b, 

where 

/Co = mis + - z'^Mdr. KLg = mis + - z Mdr, 
^ ^ 4 Jo ^ ^ ' 1 4 Jo 

ICg = ij^' Z\T)dT, 

with boundary conditions 

^(0) = Vq, z(0) = Vq, Z(0) = yg, 

Z{s) = Xq, Z{s) = Xq, Z{S) = Xg, 

and W is the Wilson loop operator 

W = (XaU^rqiVq, Xq)f^(X,)](^r,(xq, yq))i{^TM, XgU (13) 

which corresponds to the propagation of quark along the path Fg, of antiquark 
along the path F^ and of gluon along the path F^ (see Fig. ); here a, 6 are 
the colour indices in the adjoint representation. As in the case of qq system 
[17], all the dependence on the background field is contained in the Wilson 
loop operator (13). 
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The Wilson loop configuration (13) may be rewritten using the relation 
between ordered exponents along the gluon path Tg in the adjoint and fun- 
damental representations: 

^($r,(x, y))a, = (Aa)^($r,(x, y))?(A,)^($r,(y, (14) 

where the path Tg coincides with the path Tg and is directed oppositely. The 
result reads 

W = ^SpW,SpW2 - ^SpW, (15) 

where Wi, W2 and W are the Wilson loops in the fundamental representation 
along the closed contours Ci = TqTg, C2 = TqTg and C = TqTg shown at 
the Figure. 



3 Generalized area law and effective Hamil- 
tonian 

To average the Wilson loop configuration (15) over the background we use the 
cluster expansion method generalized in [18] to consider the average of more 
than one Wilson loop. For the contours Ci and C2 with the average size much 
larger than the gluonic correlation length Tg we arrive to the generalized area 
law ^ 

<>V >= ^^^exp{-aiS, + ^2)), (16) 

where a is the string tension in the fundamental representation, and Si and 
5*2 are the minimal surfaces inside the contours Ci and C2. The area law (16) 
holds for all the configurations in the qqg system, apart from the special case 
of the contours Ci and C2 embedded into the same plane, where, instead of 

(16) , one has 

<>V>=^^Y^exp(-(7(>5i->52))-(7"*>S2), Si>S2, (17) 

the string tension in the adjoint representation. The regimes (16) and 

(17) match smoothly each other at the distances between the contours Ci 
and C2 of order of correlation length Tg. 
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If the string tension is defined mainly by the contribution of second order 
correlators, then a"''^^ /a = 9/4 for the SU(3) colour group. On the other 
hand, the area law for the Wilson loop in the adjoint representation was 
observed on the lattice [19] with a^'^-'/cr ?s 2, and the same result holds true 
in the limit Nc — * oo. Having all this in mind, we assume the regime (16) to 
be valid everywhere in the qqg configuration space. 

The four-dimensional dynamics in (12) can be reduced to the three- 
dimensional one following the procedure used in [17]. Namely, choosing the 
physical time parametrization 



(r, fg), = (r, fg), = (r, fg) 



and introducing new dynamical variables 

T T ^ T ■ 

f^ii^) = 2^^°^'^^' ^^^^^ " 2i^°^'^^' ^^^^^ ^ 2S^^^^'^' 

with < T < T, we arrive to the three-dimensional representation for the 
Green function: 

G = j DfgDfgDfgDii^Dii2Dii^exp{-A), (18) 
with the effective action 

^= rrf,{|i + ^ + £i±^+^!i3+^^rrl+M+ (19) 

Jo ZfXi 2/i2 2 2 2 2 

where the surfaces 5*1 and S2 arc parametrized by the coordinates Wi^, Wi^ = 
'^i/j. ~ ~ 2- Assuming the straight-line ansatz for the minimal 

surfaces, 

Wl^ = PlZ/u + (1 - (3l)Zi^, W2n = P2Z1U + (1 - /?2)^M' 

we write out the effective Lagrangian for the qqg system as 

m2 ^2 ^^^^ ^^^2 ^^^2 ^^^2 
^"2^ + 2^+ 2 + ^ + 



7 



+(7pi £ dpi^l + ll + ap2 £ dp2\/l + ll 

- 1 ■ 

h = — Pi X {Pifg + (1 - Pi)rg), 
Pi 

h = —p2 X (/^sfq + (1 - (32)rg), 
P2 

Pl^Tq- fg, p2^fq- fg. 

To obtain the effective Hamiltonian one should define the momenta and 
express the velocities in terms of momenta. It cannot be done explicitly be- 
cause of presence of square roots in (20), and to deal with this problem the 
auxiliary field approach was suggested in [17]. However, it was shown in [17] 
that for the low values of relative orbital momenta the square roots in (20) 

— * 

can be expanded up to the second order in the angular velocities the ap- 
proximation proved to be accurate enough even for the massless constituents. 
Within this approximation the problem is reduced to the potential - like one, 
while the terms ~ can be taken into account perturbatively. The corre- 
sponding Hamiltonian in the Minkowsky space-time in the centre-of-mass 
frame is easily obtained from (20): 

2/11 2/12 2 2iip 2/xq 
Pi = P , r, P2 = -p ■ r. 

Pi + P2 Pi + P2 

where the Jacobi coordinates 

Pin + P2f2 



r = n - r2, P2 = ra 



Pi + P2 



and conjugated momenta p and Q are introduced, and /ip and /xq are the 
reduced masses 

P1P2 P2(pi + P2) 

Pp = — \ ' PQ = — \ \ — ■ 

Pi + P2 Pi + P2 + P3 

The Hamiltonian still contains the fields Pi(T), and the integration over 
{/Xj} is to be performed in the path integral representation (18) (or, equiva- 
lently, taking the extremal values of p,i in the Hamiltonian). Only after that 
the quantization should be carried out. 
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Technically, it is more convenient to proceed in a way suggested in [20]: 
first find the eigenvalues of the Hamiltonian (21) assuming fii to be c- 
numbers, and after that minimize the eigenenergies in //j. This procedure 
works with rather good accuracy for the lowest states, and reduces the prob- 
lem to the nonrelativistic three-body one, with yUj playing the role of con- 
stituent masses. We note that although the Hamiltonian (21) looks like the 
Hamiltonian of the nonrelativistic potential model, it is essentially relativis- 
tic, and the masses /li are not introduced by hand, but are calculated and 
expressed in terms of string tension a and quark masses. 

Another advantage of the above-described method is that it allows for the 
approximate solution to the problem of separating out the physical transverse 
states. Indeed, let us impose the constraint 

/i3*o - Af3(^*) = (22) 

to project out the physical hybrid state ^\ = {^q, ^), where A is the gluon 
spin index. The constraint (22) is compatible with the projector P^x (10) 
after averaging over background and introducing the variables Hi. In the 
potential-like regime one has pg = ^^fg, and we choose the physical states 
to be transverse with respect to the three-dimensional gluon momentum: 

P3* = 0, *o = 0. (23) 

We are forced to impose the condition (23), because we have neglected the 
spin dependence in the gluon Green function (9), so the condition (23) should 
be treated as variational ansatz. The rigorous analysis of the transverse and 
longitudinal gluonic degrees of freedom should be done with the inclusion of 
spin into the path integral representation for the gluon Green function. 



4 Numerical results and discussion 

In the actual calculations the Hamiltonian was supplied with the short range 
Coulomb interaction 

" 6r 2pi 2p2' ^ ' 

As a variational ansatz the Gaussian type radial wave functions were 
chosen. The constraint (23) was satisfied by taking the orbital wave functions 
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diagonal in the total angular momentum j in the gluonic subsystem, so that 
the states contain electric or magnetic gluon: 



^,— 4-4^), n-\li-^^n-UQ) + ]l^j^ymUQ)- (25) 

With this choice the electric and magnetic hybrids are degenerate, and this 
degeneracy will be removed by string corrections and by spin-dependent 
force. 

The quantum numbers of a one-gluon hybrid are given by 

p = (-1)'+^ C = (-1)'+^+! (26) 

for the states with electric gluon, and 

P^{-iy+^+\ C=(-l)'+^+^ (27) 

for the states with magnetic gluon, where / and s are the angular momentum 
and total spin in the quark-antiquark subsystem. So the possible quantum 
numbers for the ground state are 

J^c = 0^+,l^+,2^+,l^-, (28) 

where the upper/lower sign stands for the state with electric/magnetic gluon 

(26)/(27). 

The most complicated problem in the constituent approaches is not the 
relative arrangement of ground and excited states, but the absolute scale of 
masses. In the potential model large negative constant is needed to fit the qq 
spectrum, and this constant is different for the sectors with different flavour 
content. In the described approach the perimeter terms for the Wilson loop 
and/or hadronic shifts might be responsible for the constant term. We use the 
prescription that for the hybrid state with two strings the additive constant 
is twice as large as for the qq meson with only one string. 

So, the procedure used is: 1) to define the constant term for the given 
values of parameters from the fit to the S—, P— and D— wave meson levels, 
and 2) to calculate the hybrid mass with the constant multiplied by two. 

The numerical results for the spectra of hybrids with light quarks are 
listed in the Table 1 for different values of quark mass, string tension and as- 
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Slightly another procedure was used to calculate the masses of hybrids 
with heavy quarks: the constant term was taken from the fit to heavy-light 
{D— and B—) meson masses. We think that the "heavy-light" constants are 
more consistent phenomenologically for the hybrid with two "heavy-light" 
strings. The results for the ground states are given in Table 2. 

The spectra we have obtained are rather similar to the ones of the flux- 
tube model [10,11]. There is a lot of common in these two approaches, 
because both models try to account for string vibrations. Moreover, the 
Hamiltonian (21) looks quite similar to the Hamiltonian of the "one-bead" 
fiux-tube model [11]. As it was already mentioned, the numerical calcula- 
tions [11] do not support the small oscillation approximation proposed in the 
original version [10] of the flux-tube model. Hence, the constraint imposed 
that the "beads" can oscillate only in the transverse (with respect to quark- 
antiquark) direction seems to be a little suspicious. However, the heavy quark 
hybrid system was analysed [21] in the string-type regime of the Lagrangian 
(20), which matches smoothly at low j the potential-type regime described 
here, and it was shown that the effective values of api and ap2 are equal, and 
it is just the case of the "one-bead" flux-tube model. The difference comes 
from the fact that in the flux-tube the masses of constituents (including the 
bead) are fixed, while in our approach the effective masses are the variables. 
Another numerical discrepancy is due to the constant term: in the fiux-tube 
there is one string, so it is reasonable to use the constant fitted by the qq 
spectrum, while we have two distinguishable strings, and the constant de- 
fined from the qq spectrum should be multiplied by the factor of two. These 
discrepancies compensate each other, so the almost exact coincidence of the 
results seems to be to some extent accidental. 

The most important difference between the models is in quantum num- 
bers. In the presented picture the confinement is of the stochastic nature 
and is ensured by the background fields, with well-defined perturbation the- 
ory in the background. As a result, the constituent gluon carries quantum 
numbers of its own. On the contrary, the fiux-tube, being motivated by 
the strong coupling expansion, knows nothing about the gluons which pop- 
ulate the string. The excitations of the string are described by the collective 
phonon-type modes, so instead of (28) one has for the ground state 

jPC_ 0T±,1^±,2^±,1^^. (29) 

The most clear decay signature of hybrids is the suppression of a hybrid 
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decay into two S'-wave ground state mesons. This signature takes place for 
the flux-tube hybrid [22] as well as for the electric constituent hybrid [23], 
and follows from the symmetry of the wave functions involved. It means 
that P-odd hybrids (28) and (29) have the same decay properties, and the 
discrepancy between the models should reveal itself in the P even hybrid 
sector. However, all currently discussed hybrid candidates in the light quark 
sector are P odd!! 

There is growing evidence that hybrids are found at last, the belief based 
mainly on the above-mentioned signature. Indeed, the '"7r(1800) state seen 
by VES [24] decays mainly into vr/o with vrp mode suppressed; the exotic l""*" 
signal is seen in BNL [25] in the vr/i final state; the 7r(1775) is seen in charge 
exchange photoproduction [26] decaying into 7r/2 , that might be 2~+; the 
p'(1460) decays mainly into ttoi in contrast to 2 ^Si qq assignement [27]; the 
rather promising isoscalar 2~"'"?72(1870) [28] is observed with 7ra2 and ?7/o(980) 
decay models. Unfortunately, up to now no hybrid-like activity is observed 
in the P-even sector. 

5 Concluding remarks 

We have demonstrated that the perturbation theory in the nonperturbative 
confining background supports the existence of qqg bound states. The hy- 
brid mesonic excitation looks like a system of a gluon with two straight-line 
strings with quarks at the ends. For low values of gluon orbital momentum 
the problem is reduced to the potential-like one, and the resulting hybrid 
spectra are compatible with the data on light quark meson spectroscopy. 
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Table 1. Predicted masses of hybrids with hght quarks 



nig, GeV 




a, GeV^ 


Mggg, G CV 


J = 1,1^0 


i = 1,^ = 1 


i = 2,Z = 





0.3 


0.18 


1.73 


2.03 


2.2 






0.2 


1.68 


2.00 


2.18 


0.1 


0.3 


0.18 


1.71 


2.02 


2.19 






0.2 


1.68 


2.00 


2.18 





0.7 


0.18 


1.60 


1.95 


2.15 


0.1 


0.7 


0.18 


1.58 


1.93 


2.14 



Table 2. Predicted masses of hybrids with heavy quarks; 

a = O.lSGeV^, as = 0.3 



rric, GeV 


M^g,GeV 


1.2 


4.12 


1.5 


4.11 


1.7 


4.10 


rrib, GeV 


M,ig,GeV 


4.5 


10.64 


5.2 


10.64 
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Figure caption 

Wilson loop configuration corresponding to the propagation of the hybrid 
state. 
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